QUANTUM KINETIC THEORY OF TRAPPED ATOMIC GASES 



H.T.C. STOOF 

Institute for Theoretical Physics, University of Utrecht 
Princetonplein 5, 3584 CC Utrecht, The Netherlands 



Abstract. We present a general framework in which we can accurately 
describe the non-equilibrium dynamics of trapped atomic gases. This is 
achieved by deriving a single Fokker-Planck equation for the gas. In this 
way we are able to discuss not only the dynamics of an interacting gas above 
and below the critical temperature at which the gas becomes superfluid, 
but also during the phase transition itself. The last topic cannot be studied 
on the basis of the usual mean-field theory and was the main motivation 
for our work. To show, however, that the Fokker-Planck equation is not 
only of interest for recent experiments on the dynamics of Bose-Einstein 
condensation, we also indicate how it can, for instance, be applied to the 
study of the collective modes of a condensed Bose gas. 1 



1. Introduction 

The most important reason for the present interest in Bose condensed 
atomic gases, is the possibility to study in detail the dynamics of a super- 
fluid system in this case. In particular, it is possible to compare ah initio 
many-body theories for the non-equilibrium dynamics directly to experi- 
ment, which is for instance not possible for liquid helium. Two important 
issues that are of interest to us here are the dynamics of condensate forma- 
tion and the eigenfrequencies of the collective modes of the Bose condensed 
gas cloud. In both of these problems it is important to realize that the gas 
consists of two components, which, due to the harmonic confinement of 
the gas, are roughly speaking also spatially separated. More precisely, the 

1 These are the notes for a lecture delivered in July 1998 at the NATO ASI on Dynamics 
in Leiden, The Netherlands. It is a contribution to Dynamics: Models and Kinetic Methods 
for N on- Equilibrium Many-Body Systems edited by John Karkheck and to be published 
by Kluwer academic publishers b.v. 
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density profile of the gas can be viewed as a relatively narrow condensate 
peak on top of a broad thermal background. 

Theoretically, we anticipate that the noncondensate or thermal cloud 
behaves similarly as a gas above the critical temperature for Bose-Einstein 
condensation. Its dynamics is therefore accurately described by an appro- 
priate kinetic equation. However, the condensate is a macroscopic quantum 
object and the dynamics of the condensate must therefore be determined 
by an equation for its wavefunction. On the basis of these arguments we 
see that to describe the coupled dynamics of the thermal and condensate 
clouds, we need a quantum kinetic theory that is capable of simultaneously 
treating both the incoherent as well as the coherent processes taking place 
in the gas. It is the main purpose of this lecture to explain as simple as 
possible how such a quantum kinetic theory can be derived from first prin- 
ciples [1]. Moreover, as an application and illustration of the theory, we 
also present the first results on the formation of the condensate and the 
collective modes that we have recently obtained. 

It should be noted that in this lecture we use only physical arguments 
to explain and motivate the final structure of our quantum kinetic theory. 
For a detailed derivation by means of field-theoretical methods we refer to 
the existing literature [2, 3]. In addition, we mainly restrict ourselves here 
to a discussion of the so-called weak-coupling limit in which interactions 
have only a relatively small, but nevertheless crucial, effect on the behavior 
of the gas. This restriction is again made for clarity reasons only, because 
in the weak-coupling limit the theory is most transparent and therefore 
most easily understood. Moreover, once this limit is well understood, it is 
in principle straightforward to generalize and treat also the strong-coupling 
limit. However, before we can start to consider the effect of interactions, 
we first need to reformulate the theory of the ideal Bose gas in a somewhat 
unusual way, that nevertheless turns out to be very convenient for our 
purposes. 

2. Ideal Bose gas 

Let us therefore consider an ideal Bose gas in an external trapping potential 
ytrap ^ w ifa one-particle energy levels e a and corresponding wavefunctions 
Xa( x )> that can be found from the Schrodinger equation 



Using the methods of second quantization, the non-equilibrium dynamics 
of the gas is now in the Heisenberg picture fully determined by the initial 




(1) 



QUANTUM KINETIC THEORY OF TRAPPED ATOMIC GASES 



3 



density matrix p(to) and the hamiltonian 



(2) 



where (t) and (t) create and annihilate at time t a particle in the state 
Xa(x), respectively. 

If we would follow the usual treatment of an ideal Bose gas, we would 
at this point introduce the basis |{iV a }; t) in which the occupation numbers 
of all the one-particle state x»( x ) are specified, and proceed to calculate 
the probability distribution 



P({N a };t) = Tr p(t Q ) \{N a };t)({N a };t\ 



(3) 



For an ideal gas the hamiltonian commutes with the number operators 
N a (t) and the above probability distribution is in fact independent of time. 
In particular, this implies that the average occupation numbers (N a (t)) = 
Tr p(to)N a (t) are constant, which physically makes sense because without 
any interactions there is no way in which the particles can scatter from one 
state to another. 

As mentioned in the introduction, we are not only interested in the 
occupation numbers of the gas, but also in the condensate wavefunction. 
Therefore we do not want to use as a basis the eigenstates of the num- 
ber operators N a (t), but instead the eigenstates of the annihilation opera- 
tors V'a(t). More precisely, we introduce the so-called coherent states \4>;t), 
which are (properly normalized) eigenstates of the field operator 



(4) 



with complex eigenvalue </>(x) [4], and consider the corresponding probabil- 
ity distribution 



P[<f,<j>;t] =Ti p{t )\ct>-t)(<j>-t\ 



(5) 



We now need to determine the equation of motion, i.e., the appropri- 
ate Fokker-Planck equation, for this probability distribution. This can be 
achieved most easily as follows. We know that by definition the annihilation 
operators 4> a {t) obey the Heisenberg equation 



ih 



Wait) 

dt 



Tp a (t),H =e a tp a (t) 



(6) 



Since we have used the eigenstates of the annihilation operators to define 
the probability distribution P[<cf)*,(j>;i\, we also know that 

$ a (t)) = J d[<f]d[<!>] (j> a P[<l>*, 4>\ t] = (<f> a )(t) , (7) 
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where / d[0*]d[0] denotes the (functional) integral over the complex func- 
tions </>(x). Combining the last two equation we thus find that 

d 

ih—{<t> a )(t) = e a {(f> a ){t) . (8) 

Moreover, by considering the Heisenberg equation for the creation operators 
^ a {t) we also obtain 

toJ^X*) = -*«<&>(*)■ (9) 

In this manner we have thus been able to derive the equation of motion 
for the first moments of the probability distribution. However, to arrive at 
a Fokker-Planck equation for P[<j>*, (ft; t] we also need to consider the higher 
moments [5], and in our case in particular (\(ft a \ 2 ) (t) . A priori we expect this 
expectation value to be related to the average occupation numbers (N a (t)) 
and therefore that 

ih^(\cft a \ 2 )(t) = o. (10) 

Although the latter result is all that we need here, we need lateron also the 
precise relation between (\(ft a \ 2 )(t) and the average occupation numbers, 
which can be shown to be given by 

(\<f>a\ 2 )(t)=(N a (t)) + l . (11) 

A derivation of this relation is complicated by the fact that the creation and 
annihilation operators do not commute at equal times. It can, however, be 
understood physically from the fact that the second moment of P[(ft* , (ft; t] 
should contain both classical as well as quantum fluctuations. 

From Eqs. (8), (9) and (10) we conclude that the desired Fokker-Planck 
equation reads [6] 

ih^- t P[cft*, ( ft;t}= - J2JL(e a cfta)P[<ft*,<ft;t] 

+ E^Mi)*^]' ( 12 ) 

It thus contains in the right-hand side only 'streaming' terms and no 'dif- 
fusion' term. As a result there is no unique equilibrium and in fact any 
function of |</> Q | 2 is a stationary solution. Again, this makes sense physi- 
cally, because without interactions there is no way in which the occupation 
numbers can relax to an equilibrium Bose-Einstein distribution. To include 
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such relaxation into our discussion, we therefore now bring our ideal Bose 
gas into contact with a thermal reservoir. 

3. Ideal Bose gas in contact with a reservoir 

As our reservoir we take an ideal Bose gas in a box with volume V. The 
reservoir is assumed to be sufficiently large so that it can be treated in the 
thermodynamic limit. Moreover, it is in equilibrium with a temperature 
T and a chemical potential ji. The states in this box are labeled by the 
momentum Tik and equal to Xk( x ) = e' kx / 'y/V They have an energy e(k) = 
?i 2 k 2 /2m. Finally, the reservoir is thought to be in contact with the trap 
discussed above, by means of the tunnel hamiltonian 

&mt = 4r E E (*«00&(*)#(*) + £(k)&(t)$,(*)) • (13) 

Here t a (k) are complex tunneling matrix elements that for simplicity are 
assumed to be almost constant for momenta fik smaller that a cutoff Hk c but 
to vanish rapidly for momenta larger than this ultraviolet cutoff. Moreover, 
we consider here only the low-temperature regime in which the thermal de 
Broglie wavelength A = (2irh 2 /inksT) 1 / 2 of the particles obeys A » l/k c , 
since this is the most appropriate limit for realistic atomic gases. 

Due to this interaction the particles can tunnel back and forth from the 
trap to the reservoir, which results both in a shift in the energy as well as 
a finite lifetime of the state Xo( x )- The energies of the states in the trap 
therefore become complex and equal to e a + S a — iR a , where the real and 
imaginary contributions to the shift can essentially be found from second- 
order perturbation theory. Denoting the Cauchy principle value part of an 
integral by V, they obey 

and 

/r/k 
__ 5( ea + S Q -e(k))|t a (k)| 2 , (15) 

respectively. 

Introducing the retarded and advanced selfenergies HYl^ =S a ^ iR a , 
we conclude from the above that Eqs. (8) and (9) now become 

d 



fo-zrttaKt) = (e Q + " m) (MM (16) 



and 

#4 <&>(*) =-(*« + - A*) (4>* a )(t) , (17) 
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where we have also measured our energies relative to the chemical potential. 

Next, we need to consider the fluctuations, i.e., the generalization of Eq. 
(10). This turns out to be given by 

i^(|«^| 2 >(t) = -2^ a (|0 a | 2 )(t) - , (18) 

with the so-called Keldysh selfenergie equal to 

fiY% = -2iR a (1 + 2N^) (19) 

and iV^ q = l/( e ( e «+ s '«-^)/ fc sr _ j) the equilibrium Bose-Einstein distribu- 
tion function. How can this result be understood? The first term in the 
right-hand side of Eq. (18) follows simply from the fact that if we neglect 
correlations (\4> a \ 2 )(t) is equal to |(0 Q )(i)| 2 . Furthermore, the second term in 
the right-hand side of Eq. (18) guarantees that if we make use of the relation 
between (\<fi a \ 2 ){t) and the average occupation numbers N a (t) = (N a (t)), 
we recover exactly the Boltzmann equation 

§^N a (t) = -r a N a (t) + r a N? (20) 

with the correct transition rates T a = 2R a /h expected from Fermi's Golden 
Rule. 

In a similar manner as for the isolated case, we now conclude from Eqs. 
(16), (17) and (18) that the Fokker-Planck equation for our trapped Bose 
gas becomes 

ih-P[<f>*,cf>;t] = £ — -(e a + n4 + )-/i)0 a P[0*,0;t] 

d 



a U( Pa 



It is interesting to note that this Fokker-Planck equation is equivalent to 
the Langevin equation 



d_ 
dt 

where the gaussian noise obeys 



ihg- t Mt) - (ea + " /*) Mt) = Va(t) , (22) 



(r&WncV)) = l -^-h^5 a , a ,5(t - f ) (23) 
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and the fluctuation-dissipation theorem reads 

ftE* = (1 + (HEW - ftEfc")) • (24) 

Clearly, due to the fluctuation-dissipation theorem the probability dis- 
tribution P[<fi*, <p; t] for the gas in the trap now relaxes to the correct equi- 
librium state 

P[*\b °°] = n ^ + 1/2 ex P {" + 1/2 l^l 2 } • ( 25 ) 

More important for our purposes is, however, that the Fokker-Planck equa- 
tion in Eq. (21) describes simultaneously both the incoherent (kinetic) as 
well as the coherent dynamics in the gas, since it incorporates the equa- 
tions of motion for both (N a (t)) and (ip a (t)), respectively. As mentioned 
previously, this is precisely what is needed for an accurate treatment of 
non-equilibrium phenomena in Bose condensed atomic gases. 

4. Condensate formation in an interacting Bose gas 

We are now in a position to discuss an interacting Bose gas, because in an 
interacting Bose gas, the gas is roughly speaking its own thermal reservoir. 
The Fokker-Planck equation therefore turns out to be quite similar to the 
one presented in Eq. (21). To be more precise, however, we should mention 
that we aim in this section to describe the formation of a condensate in an 
interacting Bose gas under the conditions that have recently been realized 
in experiments with atomic 87 Rb [7], 7 Li [8] and 23 Na [9] gases. In these 
experiments the gas is cooled by means of evaporative cooling. Numerical 
solutions of the Boltzmann equation have shown that during evaporative 
cooling the energy distribution function is well described by an equilibrium 
distribution with time dependent temperature T(t) and chemical potential 
n(t), that is truncated at high energies due to the evaporation of the highest 
energetic atoms from the trap [10]. Moreover, in the experiments of interest 
the densities just above the critical temperature are essentially always such 
that the gas is in the weak-coupling limit, which implies in this context that 
the average interaction energy per atom is always much less than the energy 
splitting of the one-particle states in the harmonic trapping potential. 

Keeping the above remarks in mind, we find that near the critical tem- 
perature the non-equilibrium properties of the gas are in an excellent ap- 
proximation described by a nonlinear Fokker-Planck equation with time- 
dependent selfenergies fiT^ ,K (t) and effective interaction matrix elements 
V^j'd pit), f° r which in the so-called ladder approximation explicite ex- 
pressions can be derived in terms of the average occupation numbers N a (t) 
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in the gas [3]. In full detail it reads 
ih—P[<t>*,<t>;t] = 



a 0< ?<* { a', (3,(3' J 

To understand how this equation describes the formation of the conden- 
sate, we make use of the fact that in the weak-coupling limit it is appropriate 
to solve the Fokker-Planck equation with the (Hartree-Fock) ansatz 

p[4>\4>-,t] = p [ ( t>* g ,4> 9 -APiW*A'-t] , (27) 

where the complex function 4>'(x) = J2a^gXa{*-)(f) a is associated with all 
the one-particle states except the groundstate Xg( x )- Substituting the above 
ansatz in the Fokker-Planck equation, we obtain the following results. First, 
the dynamics of the noncondensed cloud is determined by the quantum 
Boltzmann equation 

^N a (t) = -TTN a (t) + r£(l + N a ) , (28) 

with a g and the scattering rates r out ' m (t) of a similar form as in the 
Uehling-Uhlenbeck equation but with a cross-section which is proportional 

t° lK!' + /3'-a (3^)\ 2 ■ Second, this kinetic equation is coupled to a Fokker-Planck 
equation for the condensate, i.e., to 

ih-P [cf>;, <p g ;t] = Qf g ( e 9 + " A* + Vg^LM 2 ) K *] 

+ ^ («P + " M + V t 9 L\^\ 2 ) PgWg, *] 

Due to the fluctuation-dissipation theorem, the probability distribution 
for the condensate relaxes to the equilibrium solution 

P [^,^;oo]«exp|-^L + ^-M+^^|^/]w 2 | • (30) 
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The Landau free energy for the condensate order parameter thus equals 



Fl\ 



(€g + Sg~ [J,)\(t)g\ 2 + 



V, 



(+) 



-to, 



(31) 



and clearly shows a spontaneous breaking of symmetry if e g +S g — \i < and 

the effective interatomic interaction is repulsive, i.e., Vg^ 9t9 > [11]. In our 
formulation of the problem these quantities are a function of time, whose 
evolution is essentially determined by the quantum Boltzmann equation in 
Eq. (28). In particular, for Bose-Einstein condensation to occur, e g + S g — \i 
has to change sign during this evolution. As a first rough calculation of 
the condensate formation, we can assume that this change of sign takes 
place instantaneously. Introducing for convenience the dimensionless time 
r = t(iY^ /8)(2Vg^g t g/kBT) 1 ^ 2 and the dimensionless condensate number 

I = \4>g\ 2 '(2Vg&g,g / IgT) 1 / 2 , this assumption leads to a typical evolution for 
the probability distribution Pq[I;t] that is shown in Fig. 1. 
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Figure 1. Typical evolution of the condensate probability distribution during 
Bose-Einstein condensation. The slowest relaxation rate is 1/to ~ 5.7. The inset shows 
the evolution of the average condensate number and its fluctuations. 
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At this point we should mention that the solution of the condensate 
Fokker-Planck equation in Eq. (29) with a fixed and negative value of 
e g + Sg— fi, is equivalent to the theory recently put forward by Gardiner and 
coworkers using different methods [12]. It even turns out to agree qualita- 
tively with experiment, although quantitatively some important discrepan- 
cies exists [13]. In our opinion these discrepancies are probably due to the 
fact that we should really solve the Fokker-Planck equation for the conden- 
sate together with the quantum Boltzmann equation for the thermal cloud. 
Moreover, and additional problem is that the experiments of interest here 
are not really in the weak-coupling limit, which substantially complicates 
the theory because more states are needed to describe the condensate. 



5. Collective modes 

As another illustration of our general nonequilibrium approach we consider 
now the collective modes of a Bose condensed gas at such high tempera- 
tures that a substantial noncondensate fraction is present in the gas. The 
experiments that have been performed at these relatively high tempera- 
tures, are in the so-called collisionless regime [14]. Physically, this implies 
that the osscilation period of the mode of interest is much shorter that the 
average time between two collisions of the atoms. Under these conditions, 
our Fokker-Planck equation for the gas gives, in the Hartree-Fock approxi- 
mation and after a transformation of Eq. (26) to coordinate space, first of 
all a collisionless Boltzmann or Vlasov-Landau equation for the long wave- 
length Wigner distribution iV(x, k; t) of the noncondensed part of the gas. 
Explicitly, it reads 

9 , T/ , x frk d , T/ , , 
— iVx,k;t + — ■ 7r W(x,k;t) 
at m ax 

" ^(^ traP ( x ) + 2 ^ (+ Mx,t))-^iV(x,k;t)=0, (32) 

where the effective interaction VW = Airfi 2 a/m can be expressed in the 
two-body s-wave scattering length a and n(x, t) is the total density of the 
gas. Furthermore, it leads to a nonlinear Schrodinger equation for the con- 
densate wavefunction (0(x))(i), i.e., to 

tftj^(x)>(f) = 

|-^ + ^ trap (x)-^ + U (+) (2n'(x,t)+n (x,t))|^(x))(t) , (33) 

with the condensate density no(x, t)) = |(0(x))(t)| 2 and the noncondensate 
density n'(x,t) = J dk AT(x, k; t)/(2vr) 3 of course adding up to the total 
density of the gas. 
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A numerical solution of the above coupled equations turns out to be 
surprisingly difficult. To nevertheless gain insight in the collective modes 
of the gas, we have therefore recently put forward a variational approach, 
in which we apply a dynamical scaling ansatz on the ideal gas results for 
both the Wigner distribution N(x, k; t) and the condensate wavefunction 
(0(x))(i). The outcome of this calculation is presented in Fig. 2, where 
we also make a comparison with experiment which turns out to be quite 
reasonable in view of the simplicity of our scaling ansatz. The main dis- 
crepancies are the two measurements halfway between the m = in and 
out-of-phase modes. This discrepancy is, however, presumably due to the 
fact that in the experiment both modes are excited simultaneously [15]. 




T/T 



Figure 2. Collisionless modes in a Bose condensed Rb gas. Curves 1 and 2 correspond 
to the m = in and out-of-phase modes, respectively. Similarly, curves 3 and 4 give 
the frequencies of the m = 2 in and out-of-phase modes. The experimental data is also 
shown. Triangles are for a m = mode and circles for a m = 2 mode. 



6. Conclusions 



In summary, we have presented a general framework in which we can dis- 
cuss various non-equilibrium problems in Bose condensed atomic gases. As 
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an example we have presented our first results on the formation of the con- 
densate and the collisionless modes at nonzero temperatures. Other topics 
that are of interest are, for example, the condensate formation in a Bose 
gas with attractive interactions, the dynamics of vortices and spin domain 
walls, the collective modes in the hydrodynamic regime, the damping of 
collective modes, and last but not least the non-equilibrium dynamics of 
trapped atomic Fermi gases. In view of the many topics that remain to be 
explored, we hope that the present lecture may motivate some of the par- 
ticipants of this summerschool to enter into this, in our view, very exciting 
area of physics. 
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